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<N : Abstract 

Neutrino-induced pion and photon production from nucleons and nuclei are important for the 
interpretation of neutrino-oscillation experiments, and these processes are potential backgrounds 
in the MiniBooNE experiment [A. A. Aquilar-Arevalo et al. (MiniBooNE Collaboration), Phys. 
^vq ■ Rev. Lett. 100, 032301 (2008)]. Pion and photon production are investigated at intermediate 

energies, where the A resonance becomes important. The Lorentz-covariant effective field theory 
contains nucleons, pions, Deltas, isoscalar scalar (<r) and vector (to) fields, and isovector vector (/?) 
fields. The lagrangian exhibits a nonlinear realization of (approximate) SU(2)l <8> SU(2)r chiral 
symmetry and incorporates vector meson dominance. Power counting for vertices and Feynman 
diagrams involving the A is explained. Because of the built-in symmetries, the vector currents are 
automatically conserved, and the axial-vector currents satisfy PCAC. The irrelevance of so-called 
off-shell A couplings and the structure of the dressed A propagator, which has a pole only in the 
spin-3/2 channel, are discussed. To calibrate the axial- vector transition current (N <-> A), pion 
production from the nucleon is used as a benchmark and compared to bubble-chamber data from 
Argonne and Brookhaven National Laboratories. At low energies, the convergence of our power- 
counting scheme is investigated, and next-to-leading-order tree-level corrections are found to be 
very small. 
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I. INTRODUCTION 



Weak pion production from nucleons and nuclei plays an important role in the interpre- 
tation of neutrino-oscillation experiments, such as MiniBooNE [l| and K2K (2)]. Pion ab- 
sorption after production will lead to events that mimic quasielastic scattering. Moreover, 
neutral current (NC) tt° and photon production produce detector signals that resemble those 
of the desired e ± signals. Finally, NC 7r° and photon production might explain the excess 
events seen at low energies in MiniBooNE. 

Ultimately, the calculations must be done on nuclei, which are the primary detector 
materials in oscillation experiments. To separate the many-body effects from the reaction 
mechanism and to calibrate the elementary amplitude, we will study charged current (CC) 
and NC pion and photon production from free nucleons. We will apply our full lagrangian 
to the many-body problem in a forthcoming paper. 

In this work, we use a recently proposed Lorentz-covariant meson-baryon effective field 
theory (EFT) that was originally motivated by the nuclear many-body problem j34l~ol]. (This 
formalism is often called quantum hadrodynamics or QHD.) This QHD EFT includes all the 
relevant symmetries of the underlying QCD; in particular, the approximate, spontaneously 
broken SU(2)l <g> SU(2)r chiral symmetry is realized nonlinearly. The motivation for this 
EFT and some calculated results are discussed in Refs. 



Here we consistently incorporate the A(1232) resonance as an explicit degree of freedom 
in this EFT, while respecting the underlying symmetries of QCD noted earlier. We are con- 
cerned with the intermediate-energy region (i^ ab < lGeV), where the resonant behavior 
of the A is important. Couplings to electroweak fields are included using the external field 
procedure [2lj], which allows us to deduce the electroweak currents. Because of the approx- 
imate symmetries contained in the lagrangian, the vector currents are explicitly conserved, 
and the axial- vector currents satisfy PCAC. Form factors are generated within the theory by 
vector meson dominance (VMD), which avoids introducing phenomenological form factors 
and makes current conservation manifest. We discuss the power counting of both vertices 
and diagrams on and off resonance and consistently keep all tree-level diagrams through 
next-to-leading order. 

The goal of this work is to use CC and NC pion production from nucleons to serve as a 
benchmark calculation. In a future paper, we will include the electroweak response of the 
nuclear many-body system to discuss pion production from nuclei. 

There have been numerous earlier studies of weak pion production off nucleons in the 
energy regime where the A is important [22 33[. It is typically assumed that the vector par t 



of the N — > A transition current is well constrained by electromagnetic interactions [29|, |31 



The uncertainty is in the axial- vector part of the current, which is determined by fitting to 



ANL [34[ and BNL [35[ bubble-chamber data. The data has large error bars, which leads 
to significant model dependence in the fitted results [3(J 32, 36[. In this work, we choose 
one recently fitted parametrization 32| and use it to determine the momentum dependence 
of the transition current vertices. We also use it to determine the constants of our VMD 
parametrization. We then compare calculations with both sets of vertices to the data at low 
and intermediate neutrino energies. 

This paper is organized as follows: in Sec. [Til we introduce our EFT lagrangian and 
calculate several current matrix elements that will be useful for the subsequent Feynman 
diagram calculations. The theory involving the A is emphasized, and the pathologies of 
introducing the A in quantum field theory are clarified, which is the basis of the lagrangian 
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construction. Then the transition current basis and form factors are discussed carefully. In 
Sec. IHIl we show the detailed calculations for CC and NC pion production and for the NC 
photon production cross sections. We initially insist on approaching this problem within 
EFT, and hence consider only the low-energy region: E^ b ^ 0.5 GeV. After that, we show 
our results in Sec. IIVI Whenever possible, we compare our results with available data and 
present our analysis. Finally, our conclusions are summarized in Sec. |V] 

In the Appendixes, we discuss isospin conventions; C, P, and T properties of various 
fields; form factors; the A propagator; and kinematics. 



II. FORMALISM 
A. Notation 

In this calculation, we use the metric g^ = diag(l, —1, —1, — 1)^, and the convention for 
the Levi-Civita symbol ^ va ^ is e 0123 = 1. The Dirac matrices are represented as (here a 1 is 
a Pauli matrix) 

and 7 5 = Z7°7 7 2 7 3 . 

Since we are going to include the A, which is the lowest N resonance, and whose isospin 
is I — 3/2, we will define the conventions for isospin indices. We will work with spherical 
vector components for the pion field, which requires some care with signs. Begin with 

A* a = T a iA A* iA , (2) 

Here a = ±3/2, ±1/2, i = ±1,0, and A = ±1/2. The upper components labeled 'a', H\ 
and 'A' furnish P^ 3//2 \ X>W, and D^ 1 / 2 ) representations of the isospin SU(2) group. We can 
immediately see that T a iA = (1, |; i, A||; a), which are CG coefficients. It is well known 
that the conjugate representation of SU(2) is equivalent to the representation itself, so we 
introduce a metric linking the two representations to raise or lower the indices a, i, and A. 
For example, A a = (A* a )* = T\ iA A iA , where iA = (§; a| 1, ~; i, A), can be written as 

A Q = T^A lA = T) B S ba 8* 5 BA A lA . (3) 

Here, 5 denotes a metric for one of the three representations. So in this convention, 
jif iA _ t* a , which is straightforward to prove. Details about the conventions are given 
in Appendix lAl 



B. Lagrangian without A(1232) 



It is widely accepted that the approximate global chiral symmetry SU{2)l <8> SU(2)r ® 
U{1)b in two-flavor QCD is spontaneously broken to SU {2) V ®U (l)s, while also being man- 
ifestly broken due to the small quark masses. To implement such broken global symmetry 
in the effective lagrangian using hadronic degrees of freedom, it was found that there exists 
a general nonlinear realization of such symmetry [37l-39|. Sometime later, in Ref. [40], the 
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concept of phenomenological lagrangians was revisited, leading to the wide use of effective 
field theory, in which chiral symmetry is realized nonlinearly. 

Here we will make use of the background field method to construct a low-energy effective 
field theory. In this method, we elevate the global symmetry SU{2) L <g> SU(2) R ® U(1)b to 
a local symmetry j2l|, 41, 



First, we will only briefly discuss how the elevated local symmetry is realized in two- 
flavor QCD, since this material can be easily found elsewhere (see Ref. (42j, for example), 
and then we show how this symmetry is realized nonlinearly in QHD. This theory is well 
developed in Refs. 0, H, [loj]. Since we take a different approach, we will detail the chiral 
symmetry realization in this section. Finally, we will talk about the electroweak interactions 
of hadrons. 



1. Chiral symmetry realization 



The charge algebra of SU(2)l <S> SU(2)r <8>U(1)b is (we will henceforth ignore the tilde 
on 5 and ?) 



ie ijk Q 



Lk 



ie ijk Qm 



[Ql , Qi] 

[Qr ) Qr] 

[Qi , Q J R ] = o , 

[Qb , Qi, R ] = , where k = +1, 0, -1 . 



(4) 



We can immediately see that massless, two-flavor QCD has this symmetry, with background 
fields including v^ 1 = v v Tj/2 (isovector vector), v£\ (isoscalar vector), a M = a* M Tj/2 (isovec- 

tor axial- vector) , s = sVj/2 (isovector scalar), and p = pVj/2 (isovector pseudoscalar), 
where i = x, y, z or +1, 0, —1: 



£> = £qcd + + Bv^ s) + 7 5 a M )g - q(s - ij 5 ?)q 

= C QCD + IlI^ + B^ s) )q L + q Rl ^ + B^ s) )q R 
~ 9l(s - ip)qR ~ q R (s + ip)q L ■ 



(5) 



- a**, q L 



Here, r M = v M + a M , 
baryon number. The symmetry transformation rules are 

B 



|(1 - 7 5 ) q, qR = |(1 + 7 5 ) 9, and 5 = 1/3 is the 



qiA 


-> 


exp 


.9{x)~ 
3 


^exp 


qR 


-> 


exp 


.0(a;)~ 
3 


exp 




-)■ 


L PL* + iL d^L ] , 




— >■ 




+ iRd»tf , 


I 1 - 

V W 


— >■ 


u 






s + ip 


-)■ 


R(s + ip)L^ , 


s 



iOrAx) — 



q LB = exp 



-iO Ri (x) 



q R = exp 



.0(x) 



,fl(z) 
3 

RqR 



s — ip — > L(s — ip)R' ( 



(6) 
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We can also construct field strength tensors that transform homogeneously: 



jhysn = d^ljj — d v l^ — i [1^ , l v ] — > hf^l) , 



(7) 



Meanwhile, to conserve C, P, and T symmetry, we have the corresponding transformation 
rules shown in Appendix iBl 

Now we proceed to discuss low-energy nuclear theory involving 7r*, p 1 , N A , and the chiral 
singlets and <p Q. As noted earlier, chiral symmetry is spontaneously broken at low 
energy in the chiral limit, and the symmetry is realized nonlinear ly: 



U 



cxp 



i = Vu 



U = exp 



.7T< 

'a 



hdfjv 



[IV 
[IV 



(dp + iu M - iv( a ) M S)^Vi9 ->■ ex P H0(aO-B] hi(d^) 



(8) 



In the preceding equations, t l are the generators of reducible representations of SU{2). 



Specifically, they could be generators of V^ 2 ^ (BV^ '(BT> [ £' z> , which operate on non-Goldstone 
isospin multiplets including the nucleon, p meson, and A. We will generically label these 
fields by ip a = (Na, Pi, A fl ) a . Most of the time, the choice of f is clear from the context. 
B is the baryon number of the particle. The transformations of the isospin and chiral 
singlets and <fi are trivial. We will also make use of the dual field tensors, for example, 
p ± iiv ^ e nva/3 ^ w hi c i 1 have the same chiral transformations as the ordinary field tensors. 
Here we do not include the background fields s and p mentioned in Eq. (JSJ), which are the 
source of manifest chiral-symmetry breaking in the Standard Model. 
The realizations of C, P, and T symmetries are given in Appendix O 



,(3/2) 



2. Power counting and the lagrangian 

Based on the transformation rules of the building blocks listed above, we can begin to 
construct the low-energy EFT lagrangian. The organization of interaction terms in this 
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lagrangian is based on power counting [5j. The power counting in EFT essentially assumes 
that when the interaction structure becomes more complicated, i.e., more fields and more 
derivatives are introduced, its contribution to physical observables becomes less important. 
Similarly, loop contributions will be suppressed more when the number of loops gets big ger . 
The validity of this assumption is connected with Naive Dimensional Analysis (NDA) [43J, 



44j | which assumes the strength (coupling) of the interaction is of order unity when the 
appropriate dimensional scale factors have been included. This "naturalness" can be checked 
only after the calculations are finished. 

To make the power counting transparent, we can associate with each interaction term an 
index 

u = d + - + b. (9) 

Here d is the number of derivatives (small momentum transfer) in the interaction, n is the 
number of fermion fields, and b is the number of heavy meson fields. Since the lagrangian is 



well developed in Refs. [10|, |45[, we just outline the lagrangian here. We begin with 
Cn(p&) = N(i'f[dn + l 9 P Pn + WvVp] + gAl^a^ -M + g s 4>)N 

+ ^-NN Tr(aX) + ^ Tr (W) 

+ ^Na^(2X^f s , u + A«F«)iV , (10) 
~ ~ «- 

where is defined in Eq. (JHJ), d v = d u — (d u — iv v + iv( s )„), and the new field tensors are 
= d^Vv - and 

Pfiu = d[^p u] + ig p [p^ , p v ] + , p v ] - p, -B- v) ->• h p^ u h ] . (11) 

The superscripts and ^ denote the isospin. 
Next is a purely mesonic term: 

£ m eson(^4) = i«9^^0+^^Tr[^(^f/)t] + i/>2 Tr(f/ + f/ t_ 2) 



+ 1 + n„ m n it (Pupr) — — — r — r — =- m:0 

\ 'P A/T P \rPr I n 1 O J/ /I /I /T2 / 



M y p w ' V 2 3! M 4! M 2 

+ A- + \ ffV^ . (12) 

The v = 3 and ^ = 4 terms in £ppn(i><4) are important for describing the bulk properties 
of nuclear many-body systems p, |4y, |47j . The only manifest chiral- symmetry breaking is 
through the nonzero pion mass. 
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Finally, we have 



iV,7r(i>=4) " 

- Wl Tr(F« O + w 2 Tr^F^) + u; 3 Tr (a„z [a u , ) 
-<W^Tr(p^) 



S/ll/ 



+ C j^fN^NTr (W~p, v ) + e -f^Nl^NV, u . (13) 

Note that Cn,tt(u=a) is not a complete list of all possible u — 4 interaction terms. However, 
/?(°) and will be used in the form factors of the nucleon's vector current, 0^2,3 will 
contribute to the form factor of the pion's vector current, and g p7T7T will be used in the form 
factors that incorporate vector meson dominance. The constants Ci,ei,Ci p , and e\ p will be 
explained later when we discuss photon production. 

3. Contributions to current matrix elements from irreducible diagrams 

By comparing Eq. (jSJ) with the electroweak interactions of quarks in the Standard Model 
we can determine the form of the background fields in terms of the vector bosons 



W±, Z„, and A^. 



r° g t° 

lu = — e — A u H — sin 2 6 W — Z u 

M 2 M cos9 w 2 M 



-e — A, H — sin 2 6L — Z„ , (15) 

2 cosS,,, 2 v ; 



Vf s i„ = —e-A u -\ ^— sin 2 6L - Z u , (16) 

(sJ/i 2 cos^ 2 M v ; 

where (7 is the SU(2) charge, and 6 W is the weak mixing angle. Furthermore: 
r° g r° r° 

= " e T Am + ^ sin2 ^ ~2 Zm ~ co^T 

+ interference terms including (WZ) , (WA) , (WW) , but no (Z.A) , (17) 

r° g r° 

/ii^ = -e — A M + ^ - sin 2 6 W — Z M (no interference terms), (18) 

Uu = ~e ~ A M + sin 2 9 W - Z M . (19) 
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If we define [see Eq. ((Sj)] 



-eJ* M A» 



-^-J^ C Z» - gV ud J L +l ,W +l » - gV ud J^W^ 

COS 1/7;; 



(20) 
(21) 



and use Eqs. f lT4|) to ffTB"]) . we can easily discover 



4 = nO^ + A 



H/lJ j 



— 



2 \Viii Aifi) , 



J, 



v.. 



- J, 



B 



J, 



J M - sin ^ 



(22) 

(23) 

(24) 
(25) 



Here, is the baryon current, defined to be coupled to v^. These relations are consistent 
with the charge algebra Q = T° + B/2. (B is the baryon number.) V 1 ^ and A lfl are the 
isovector vector current and the isovector axial- vector current, respectively. We do not 
discuss "sea gull " terms of higher order in the couplings because they do not enter in our 



calculations Il0j,l50j 



Based on the equations given above and the EFT lagrangian, we can calculate the matrix 
elements (N\V*, A^, Jj?\N) and ( N ; 7r|Vj, A^, J*\N) at tree level; loops are not included. 1 
Since non-Goldstone vector bosons are included here, then by vector meson dominance 
(VMD), we can extrapolate the current away from Q 2 = to some extent [l(J 20]. The 
results are given below, and the explicit calculations are shown in Appendix |D] Note that 
q^ is defined as the incoming momentum transfer at the vertex; in terms of initial and final 



nucleon momenta, q^ = p^ — p^ 



Similarly, q^ + p, 



Pnf + ^7r f° r pi° n production. 



(N,B\V'\N,A) 



/ n \ T "\A\— ( 1 ^^T?V,md ( ^l^~ 4% . nI pV,md 0~v, v iq" 

(B\ — \A) u f ( 7^ + 2bF l h 2F 2 — — I //, 



2M 



(B\ — \A) u f r Vfl (q)ui 



° n i cU (B\ T -\A) u f ( 2F?>™% + 2F™ 



2M 



(26) 
(27) 

(28) 



1 The expressions for the currents listed below differ from those in Refs. 10, 5l| because contributions from 
non-minimal and vector meson dominance terms are included here. 
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{N,B\tf\N,A) = 5£u f ^ + 25F{ 



s,md q 2 lfi- 4% _|_ 2pS,md cr f iuiq L 



f ' -' - 2M 1 " ; 



on shell 



It; 



(29) 
(30) 
(31) 



(N,B;n,j,K\A;\N,A) = --f {B\ T -\A) urf Ui 



x 



• \ '.mil, ,_ ^2^5 , ' \Q ~ k~ ).'////' ' ('I ' K ) /i'Ii- 



g llu + 25F 1 v ' ma ((q-k 7T ) 



2A"> 



e l . k T k 

-J- (B\-z-\A) u f T AnfM (q, k w )ui . 

Jtt * 



Here (m p = 0.776 GeV, m v = 0.783 GeV): 



V.md 



-,V,md 



U 1 + ^ q ,_9_ R _ 



2 V g~, q 2 -m 2 



, = -1-35, & = 2. 
) , A« = 1.85, f p = 3.04 , 



S.md 



^S.md 



Id + aHlj-til q2 ) AW = -140 ^ = 395 

^^V) , \ (0) = -0.06, / w = -0.19 , 
3# 7 q 1 -ml) 



2A<°> 



SF V/S,md^ = F V/S,md^ _ pV/S.md^ 



5F 



V/S,md ( 2 



f S ' md {q 2 ) - F? /S ' md (0) . 



(32) 
(33) 



(34) 
(35) 
(36) 

(37) 

(38) 
(39) 



If we follow a procedure similar to that used in calculating (iV|V^|iV) and (iV; 7r|A^|iV), 
we can expand the axial-vector current in powers of q 2 using the lagrangian constants qa 
and (3 A l \ In fact, we can improve on this by including the axial- vector meson (a llx ) contri- 
bution to the matrix elements, which would arise from the interactions: g ai N^rfai^N and 
c ai Tr (F^^ai^) . Here = clu^t 1 /2 and = d^a lu — d u ai^, where are the fields 
of the ai meson (whose mass is denoted as m ai = 1.26 GeV). Then we obtain (details can 



9 



be found in Appendix ID]) 

(N,B\A^\N,A) : 



T 



<W) (B\-\A)u f [^ 



% 4 



mi 



, T 



(B\-\A)u f Y A M^ 



(40) 



(N,B,*,j\V*\N,A) = ^(B\-\A)u f [GT(0)w 



-mid i 



fn ^ ' 2 

+ 5GT((q-K) 2 



q-(q- K)g^u -(q- K)^ 



(q - k n f 



Yi 5 \u 



6 k T 
= -j- ( B \-^\ A )ufT Vntl (q, K)m 

Here, the definitions of G^ d (q 2 ) and 5G A ld (q 2 ) are 



(41) 
(42) 



r md (n 2\ _ „ <? 2 2c ai 5f ai g 2 



1.26, 



9a 

For the pion's vector current form factor 

(vr, fc, MKK, j, K-q) = ie\(2kv - q)^ ( 1 



M 2 — 

2.27, c ai (? ai = 3.85 



G7V) - GT'(O) = G^V) 



9pinT q 



9a 



(43) 
(44) 



g 1 q 2 - m 2 



pion on shell 



ie\ (2K ~ q)»F™\q 2 ) 
(pion off shell) 



QpTVTl 

9y 



1.20 



(45) 



(7r,k,k w \V^7r,j,k n -q) = ie 



XI 



h — 



q ■ fc„ 



5m 



with 8F™ d {q 2 ^ 



it^Pvuiqih) , 
F™V) - F™ d (0) . 



(46) 
(47) 



To determine the couplings in Eqs. ( 134")) . (!35|) . (|36|) . ( 137|) . (143]) and (l4"5l) . we compare 
our results with the conventional experimentally fitted form factors [5j, |52j. We make the 
behavior of the form factors near q 2 = as close to the conventional form factors as possible. 
In this calculation, we fit our 'md' form factors to those in 52] for the nucleon's vector and 
baryon current. The conventional nucleon's axial- vector current used to fit our G 



md 



IS 



parameterized the in literature [53| as Ga(q' 2 ) = 9a/ (1 ~ q 2 /M A r ) 2 , with qa = 1-26 and 
Ma = 1.05 GeV. As shown in Ref. 20), the form factors due to vector meson dominance 
become inadequate at Q 2 ~ 0.3 GeV 2 . This is also true of the axial parametrization. This 
indicates that the EFT lagrangian is only applicable for Ei ^ 0.5 GeV in lepton-nucleon 
interactions, above which Q 2 exceeds the limit. This will be clarified in the kinematical 
analysis of Sec. IIV Al 
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C. Lagrangian involving A(1232) 



1. Chiral symmetry and power counting 

The A* a belong to an I = 3/2 multiplet as non-Goldstone particles in the low-energy the- 
ory. The chiral symmetry realization involving non-Goldstone particles, ip a = (Na, Pi, A a ) a , 
has been given generally in Sec. Ill B 11 Moreover, in the power counting of vertices, the A 
is counted the same way as the nucleon. 



2. Spin-3/2 particles as manifest degrees of freedom in field theory }54j 



We briefly discuss the well-known pathologies of high-spin particles in field theory [55M59 
It was discovered that with strong couplings, or strong fields, or large field variations, a field 
theory involving high-spin fields cannot be self-consistent. These investigations were carried 
out in the lagrangian formalism with a finite number of interaction terms. The reason for 
the pathology is that the unphysical degrees of freedom of the A may contribute when the 
constraints on the A fields are lost after adding interactions. However, it has more recently 
been realized that the number of degrees of freedom is correct in low-energy effective field 
theory, which has, in principle, an infinite number of interaction terms, as long as we work 



in the limit of low-energy and weak (boson) fields |6CH62| . As we know, in the low-energy 
effective field theory, the interactions and boson fields are scaled by 1/M « 1/(1 GeV). The 
essence of the argument is that in the perturbative picture, the spin- 1/2 components of 
the off-shell A (treating the A as a stable particle) will behave as local contact interactions 
without pole structure, which should be calibrated together with the complete set of contact 
terms in the effective lagrangian. 

Another issue is about the so-called off-shell couplings, which have the form 7 / /0'\ 
d^ip^, ip Jfi, and d^ip , and which had also been discussed together with pathologies. From 
the modern effective field theory viewpoint, it has been concluded 63, that these off-shell 
couplings are redundant. The physical picture of this will be clarified after introducing the 
lagrangian in Sec. Ill G 31 



3. Lagrangian and A propagator renormalization 

Consider first £a;7t,p,v,</> ip ^ 3), which is essentially a copy of the corresponding la- 
grangian for nucleons: 



L^v* = y ,(i@-h P i6-h v y-m+ h s <j>)} b a A bv + h A ^ ^ 7 5 A£ 

- ^A x p^A' - 4^ A A V^A* 
4m Am 

- ^A x v^A x + ^ A A A A Tr (a^a M ) . (48) 
m m 

Here the sub- and superscripts a,b = (±3/2, ±1/2), so the isovector mesons use the isospin- 
3/2 representation. A few words on pathologies and off-shell couplings are in order. Based 
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on CA-,TT, P ,v,<j>, we can find the free propagator for the A. Details are shown in Appendix [0, 
and here we ignore \w indices and isospin indices: 



-i 



-.O(l-L). 



p(§). 



m 



+ pd 1 -) 



le 
1 



m 



_x 12 



rn 



Pd x ) . (49) 



The operator projects the general Rarita-Scwhinger field tp^ to spin-3/2 objects, while 
P^- 1 ) = 1 — pb) is the orthogonal projection operator. In Appendix [Ej we also show that 
the self-energy £^ = £ A g MJ , + ST,^ can be written as 



p(i)£ A p(!) + p(f^)sp(f ± ) 



(50) 



We can then renormalize the A propagator as follows: 



S, 



+ (^ + S J(f-L))( E (§) +E (i±)) (5 - 



0(f) 



+ 



+ 5 



O(f-L) 



+ 5 



o(f±) s (|±)^o(|±) 



+ ■ 



(51) 



From this we can conclude that the renormalized propagator Sp 



S { } ] + S^. The 



(I) 



resonant contribution is S F 2 



5 o(§) + ^(DwSlcd) 



Q 2 



(f-L) oO(l-L) 



Sp 2 



Y*(^S F 2 '. The background contribution is 



+ S F K * W^Sp 2 J . We can see that renormalization will shift the pole 
position of the resonant part. Moreover, as long as power counting is valid, i.e., 0(S/M) <C 
1, we will always be far away from the unphysical pole in the renormalized non-resonant part 
proportional to {1/[1 — 0(S/M)]} _1 . This also suggests that we will not see an unphysical 
pole in the renormalized propagator, when we work in the low-energy perturbative region, 
where power counting makes sense. So we indeed have the right degrees of freedom in our 
low-energy theory, namely, a single pole at the resonance. Thus perturbative unitarity is 
not obviously violated in this theory with high-spin fields. 

Another issue is the 1/p 2 singularity in the projection operators. In principle, when we 
are in the low-energy region, p 2 = (p n i + k) 2 is always positive (timelike) in all the channels. 
Here k is some general small momentum compared to p n i, which is the nucleon's momentum 
that is almost on shell. 

The preceding discussion also helps to clarify the redundancy of the off-shell couplings. 

(-) 

The self-energy due to these couplings will not contribute in the renormalizaton of S F 2 , but 
it will indeed change the non-resonant part. However, the effects due to these couplings 
can be expanded in powers of the momenta. So they will essentially look like higher- order 
contact terms without the A. This justifies the redundancy of these couplings. To ignore 
these couplings in a way that does not break term-by-term chiral symmetry, we can always 
combine <9 M with pion fields so that it becomes d^. This indicates that those couplings with 
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or 7^ contracted with A M can be ignored without breaking manifest chiral symmetry 

For a concrete example of these results, see Appendix [0 
To produce the N ■(->• A transition currents, we construct the following lagrangians [y ^ 

4): 

C A , N;7r = h A A afl T^ A a lf ,N A + C.C. , (52) 



.4 



A^T^^F^)^ + C.C. , (53) 



M 2 

Ca;N; P = ^A; %7 5 Tt^^ A + ^A;, 7 5 Tt^(^),iV A 

+ S A V pH^a + C.C. (54) 

Terms omitted from these lagrangians are either redundant or are not relevant to our cal- 
culations |50(. 



4- Transition currents 

It is easy to expect the validity of the following definitions: 

(A,a,p A \V^\N,A,p N ) = T^ A u Aa (p A )T^(q)u N (p N ) , (55) 
(A,a,p A \A^\N,A,p N ) ee Ti iA M AQ ( PA )r7(g)^( Pw ) . (56) 
Based on the lagrangians given previously, we find (note that cr liu e fiual3 oc 2<r a/3 7 5 ) 

rr = (g a r- ^h 5 + 2 -^0 1 (?V - jTsV 



8c 6A (g 
M 2 



2\ 



(f^'q.-f , (57) 



FT = -h A (> - + ^ ( 5 V - 9 «<V ) - ( 5 V - 9 



O/i 



- ^ 9 Wg„ , (58) 

Q A (g 2 ) = CiA + — 2 J' « = 1,3,6, (59) 

2# 7 q 2 -m 2 

c 1A = 1.21, c 3A = -0.61, c 6A = -0.078, 

^£ = -4.58, ^ = 2.32, ^ = 0.30. (60) 

9-y 9"/ 9"/ 
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Quite similar to the qa(q' 2 ), we can also introduce axial- vector meson exchange into the 
axial transition current, which leads to a structure for the di^,(q 2 ) that is similar to the 
vector transition current form factors. There is one subtlety associated with the realization 
of liA(q 2 ), which is the same as the one detailed in Appendix [D] for G™ d (q 2 ): with our 
lagrangian, we have a pion-pole contribution associated only with the h A coupling, and all 
the higher-order terms contained in Sh^q 2 ) = h^iq 2 ) — conserve the axial transition 
current. With the limited information about manifest chiral-symmetry breaking, we will 
ignore this subtlety and still use the form similar to the c\&(q 2 ) to parameterize liA{q 2 )- The 
axial-vector meson couplings h& ai and diA ai will be combinations of g ai and the coupling 
strength of Aa\N interactions. So we have 



h A {q 2 ) 

diA(q 2 ) 
h A 



9 9 ' 

Q ~ m ai 



diA + diAax 



q 2 - m 2 ai 



1.40, 



d 



2A 



-0.087, 



2,4, 7, 



diA = 0.20, G?7A 



-0.04, 



-3.98, d 2Aai = 0.25, d iAai = -0.58, d 7Aai = 0.12. 



(61) 
(62) 

(63) 



To determine the coefficients in the transition form factors shown in Eqs. fl60l) and fl63l) . 
we will compare ours with the conventional ones used in the literature. In Refs. 30, |32| for 
example, the definition is 



(A,-|i^ + |iV,-^) = u a (pa ) 



+ % (q ■ pn g a " - g a P%) 



7 



+ 



5 9 M 2 



u{p N ) 



(64) 



u{p 



Ni 



(65) 



where Yy 1 and Y^ 1 are defined in Eqs. fl55l) and fl56|) . The basis given above is known to be 
complete for the transition matrix element. The phenomenological form factors are listed 
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below 32 



V/2\ 



Here 



a 



G D {q 2 ) 

C?(q 2 ) 
C£{q 2 ) 

C£{q 2 ) 



2.13 



1 - (g 2 /4M 2 ) 

-1.51 

1 - (g 2 /4M 2 ) 

0.48 

1 - (g 2 /0.776M 2 
1 



[1 - (q 2 /M 2 )} 



2\1 2 ' 



G D (q 2 ) , 
G D {q 2 ) , 

G D (q 2 ) . 

and M v = 0.84 GeV. 







4 

C£{q 2 



C£{q 2 ) 



M 2 



m 2 — q 2 



1.14 1 + 



1.21g 2 



2GeV 2 -g 2 ; [1 - (g 2 /M 2 )] 2 



(66) 
(67) 
(68) 

(69) 

(70) 
(71) 

(72) 

where M A = 1.29 GeV. 

(73) 



To equate the two different representations of the transition currents when q 2 = and 
the baryons are on shell, we have: 



ClA 
C3A 
C6A 

h A 

d-2A 
G?4A 



Cl m-M(CY + CY) 
2 2M 2 

3 (CY - CV) 

2 4 

3 (CX + CD 
2 16 

3 



a 



5 ) 



3 qt_ 

2 4 ' 
3 , m + MCt 



d 



7A 



2 V 2 

3 C£ 
2 8 



+ 



2M 2 



(74) 

(75) 
(76) 
(77) 
(78) 
(79) 
(80) 



To determine our own form factors, we assume that the relations above hold not only 
when q 2 = 0, but also in kinematic regimes with finite q 2 . It can be shown that at low energy, 
the differences in observables due to using the two bases, with these relations applied, are 
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negligible. This shows that our basis of invariants is also complete and partially justifies 
extending the preceding linear relations to kinematics with finite q 2 . Alternatively, all of the 
q 2 dependence of these qa and di& form factors can be realized in terms of meson dominance. 
We then require that the meson dominance form factors be as close as possible to the ones 
produced by the conventional form factors in Eqs. (1741) to (1801) . However, when we compare 
our meson dominance form factors with those in literature, we clearly see the inadequacy of 
the leading-order meson dominance expressions above Q 2 ~ 0.3 GeV 2 . 

III. FEYNMAN DIAGRAMS 




(a) (b) (c) (d) (e) (f) 

FIG. 1: Feynman diagrams for pion production. Here, C stands for various types of currents 
including vector, axial-vector, and baryon currents. Some diagrams may be zero for some specific 
type of current. For example, diagrams (a) and (b) will not contribute for the (isoscalar) baryon 
current. Diagram (e) will be zero for the axial- vector current. The pion-pole contributions are 
included in the vertex functions of the currents. 

Tree- level Feynman diagrams for pion production due to the vector current, axial- vector 
current, and baryon current are shown in Fig. [U In this section, we will begin to calculate 
different matrix elements for pion production and photon production. The Feynman dia- 
grams for photon production can be viewed as diagrams in Fig. [1] with an outgoing n line 
changed to a 7 line. It is easy to see that diagram (e) in Fig. [I] will not contribute to photon 
production, since there is no vertex connecting a pion and a photon. 
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A. Renormalized A propagator 



Since the one-pion-loop self-energy has already been calculated in Ref. 45], we simply 
quote the result here: 



H (y + m)Pi| ^ + 0(r/m) x non-pole terms, (81) 

3m 2 

F « = 12^(sfp (p2 + M2 + 2Mm) 



X 



[(p 2 - M 2 ) 2 - (p 2 + 3M 2 )m 2 ] a/ (jo 2 — M 2 ) 2 — 4p 2 m 2 . (82) 



However, to make the calculation simpler, we will set II = 0. We take m = 1232 MeV as the 
Breit-Wigner mass [65]. Note that T is implicitly associated with a factor of Q[p 2 — (M + 
m n ) 2 }. 

A few words on the 1/p 2 singularities in the projection operators are in order here. (See 
Appendix [El) When the A is in the s channel, we never see these singularities, since p 2 is 
timelike. When the A is in other channels, forp 2 ^ (m n + M) 2 , we never see the singularities. 
For p 2 ^ (m n + M) 2 and hence T(p 2 ) = 0, the apparent singularities are actually canceled 
out in our approximation scheme. It can be easily checked that if we set II = 0, then when 
T = 0, Sp U (p) — > Sp^ip), so no singularities will appear. This is the advantage of setting 
II - 0. 

We have to remember, however, that the whole calculation is valid only in the low-energy 
limit, and in this limit, we would not see a A in the u channel that is far off shell, because 
both q and k n are tiny. Hence the 1/p 2 singularities should not be a problem in low-energy 
effective theory from a very general perspective. 



B. Power counting of Feynman diagrams (v) [6 

First let's outline the calculation of the interaction amplitude M. Consider CC pion 
production (in the one- weak-boson-exchange approximation): 

M = AV2G F V ud (J^)(J^ ad)itl ), , (83) 
where i = +1,-1. In Eq. ([8*31 . Gf is the Fermi constant, V u d is the CKM ma- 
trix element corresponding to u and d quark mixing, (Jp^) = (l(J)\ JunWii^i)} ^ an d 
^j(had)in^ ^ (NB, irj\ J]f| NA). (Jua ) * s wen known, so in the following, we focus on 
calculating ( Jr' iad - >lM } 7r . 



L 

For NC pion production: 



M = AV2G F (4 e Si)(J { N^ , (84) 



where (J^cl) is the well-known leptonic neutral current matrix element, and (J^c }n = 
{NB,7rj\J^ c \NA). For NC photon production, we have a similar expression: 

M = AV2G F (4 e Sl)(J ( NC^ , (85) 
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where (J^) 7 = (NB n \J* c \NA). 

Now consider the power counting for ( j( had ^) 7r ( 7 ) in Eqs. flBBl . and fl85|) . The order 
of the diagram (i/) is given by normal power counting jisl j: z/ = 2L + 2 — 1 + — 2), 

where L is the number of loops, i? n is the number of external baryon lines, V{ = di + | rii + 6; 
is the order of the vertex (z>) defined in Sec. Ill B 21 and #j is the number of times that 
particular vertex appears. 

As pointed out in Ref . [66| , by including A resonances in calculations, we have a new mass 
scale 5 = m — M « 300 MeV. We must also consider the order of the A width V. Formally, 
it is counted as 0(Q 3 /M 2 ); however, numerical calculations with Eq. ( 182]) indicate that it 
should be counted as 0(Q 3 /M 2 x 10). Because of these two issues, we have to rethink the 
power counting of diagrams involving 5 in two energy regimes. One is near the resonance, 
while the other is at lower energies, away from the resonance. In the resonance region, the 
A propagator scales like 

F iT \MJ 10t 0{Q 3 /M 2 ) iO{Q 2 /M) 0{Q)iO{Q)' 1 ; 

where the 0(1/M) comes from non-pole terms. In the lower-energy region, 



2[5 - 0(Q)} - 10z 0(Q 3 /M 2 ) \M J 0(Q) 5 \M J 0(Q) M 

(87) 

So compared to the normal power counting mentioned above, in which the nucleon propa- 
gator scales as 1/0(Q), for diagrams involving one A in the s channel, we take v — > v — 1 
in the resonance regime and v — > i/+ 1 away from the resonance. This partially justifies the 
strategy of incorporating non-resonant diagrams at low energ ies, while ignoring them in the 



resonance region when fitting the form factors in this region [3JJ, |32 



C. CVC and PCAC 



We will calculate matrix elements of currents and test the conservation of the vector 
current and the baryon current, and also partial conservation of the axial- vector current. 



1. Diagram (a) 

Diagram (a) in Fig. [TJ leads to a vector current (k n is the outgoing pion's momentum) 

(v*% = ~ t% T\ iA u f k* s FXa (p) iT(p; g.ftte • (as) 

Jit 

Here Yy i {p- 1 q,Pi) is defined in Eq. (j551) . Momentum conservation gives p = q + Pi, and 
T a jB Tl iA = 5)5* - \{t 3 t*)£ = 18)8* - |ef {r k ) A B = \5)5 A + \e\ k (r k ) A , where the subscript 
j denotes the isospin of the outgoing pion. Vector current conservation is obvious, and 
v nr ^ 3 in the lower-energy region, while v r ^ 1 in the resonance region. Here the higher- 
order terms in v come from including form factors at the vertices. 
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The axial- vector current matrix element is 



<^*% = --^T a Bj T^ A u f k x w S FXa (p)TT(p;q,p l )u l . (89) 

JIT 



Here r^ M (p;g,p;) is defined in Eq. ( 156|) . PCAC is also obvious, if we check the structure of 
T^ 4 , and v nT ^ 2, i/ r ^ 0. 

The baryon current matrix element is 

(JS)* = • (90) 

Now we examine the NC matrix element (^jv^ } 7 - First, based on the relations given 
in Eq. f[2"o"j) . we define 

FtffaQiPi) = (^-sin 2 ^)r^(p;g,p i ) + ir^(p;g,p i ) , (91) 

r N a (p f] q,p) ee j°T^(p; -? lP/ )7° , (92) 

(A,a,p|J^ c (g)|iV,A^> = Ti 0A w Q (p)r^(p;g, ft )«fe) , (93) 

(JV.AP/l^o^l^^P) = T a 0A u(p f )T^(p f ;q,p)u a (p) . (94) 
Then we find [k is the outgoing photon's momentum and e* x (k) is its polarization] 

(r NC ), = eT^rf^^e^^r^^j-A;^)^^)^^;^^)^ • (95) 



CVC and PCAC are straightforward to verify here. For the vector current, v nr ^ 4, v r ^ 2, 



while for the axial- vector current, v nr ^3, z/ r ^ 1. 



2. Diagram (b) 

Diagram (b) in Fig. [TJ leads to the vector current 

(V*% = -^nTl A u f r v a (p f ,q,p)S FaX (p)ktu l . (96) 

Jtt 

Here Ty a (p f ;q,p) = 7°F{^(p; -q,Pfh°, P = -q+Pf, and T%T\ A = 5]5 A - KtV^b = 
| ~ t e% j\Sj k ^B- The conservation of the vector current is obvious, and u nr ^ 3. 



The axial-vector current matrix element is 

= -^T-Tl A u f r A a (p f] q,p)S FaX (p)ktu t . (97) 

Jtt 



Here T A a (pf] q,p) = 7°r^ t (p; —q,Pf)j°, p = —q +Pf, PCAC is again obvious, and v nr ^ 2. 

The baryon current matrix element is zero ((J^) n = 0), and the NC current matrix 
element for photon production is 

( J nc)i = eT^T t a A u f T f ^{pf,q,p)S Fa p{p)T^ X {p-,-k,p i )e* x {k)ui . (98) 

Both CVC and PCAC are obvious here. For the vector current, v nr ^ 4, while for the 
axial- vector current, v nr ^ 3. 
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3. Diagrams (c) and (d) 

These two diagrams lead to a vector current 



u {B\ T -f-\A)u } ^S F {p)^ v { q ) Ul 



T T. 



- ^ {B\^\A)u f T» v { q )S F {p) fin! m 



(99) 



Here S F (p) = S F (q + pi) is the nucleon propagator, r^(g) has been defined in Eq. fl27|) . and 
v ^ 1. To prove CVC, one must consider diagrams (c), (d), (e), and (f) together. 
For the axial- vector current, we find 



Wa 



(B\^\A)u f T^q)S F (p) 



2 2 



(100) 



Here, T^(g) has been defined in Eq. ( 1401) . PCAC is obvious, and v ^ 1. 
For the baryon current we have 



l jA(Bq\A)u f M 5 S F (p)TMu t 



Wa 



-f(B\j\A)u f T%( q )S F (p) %r 



U; 



(101) 



Here, T^(g) has been defined in Eq. f l30j) . It is easy to see that the baryon current is 
conserved and that v ^ 1. 

Finally, for NC photon production, we get 



'Ah = en/eK*) ((y)B C r^(-A;) + ^r^(-^^(p) 



r \ A 

~2 )c 



;i-sin 2 W( g ) + lrM 



5 A 



c 

IB 



;i-sin 2 w( g ) + irM 



^sin 2 ^r£(g)) 



xS F (p)el(k) 
where we use the shorthand 



6, 



: y )b = (b\ y \a) . 



(102) 



(103) 



One can verify the conservation of the vector current and the baryon current, as well as the 
partial conservation of the axial- vector current. For all three currents, power counting gives 
2/^1. However, this naive power counting does not give an accurate comparison between 
the A contributions and the N contributions at low energies, as we discuss below. 
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4- Diagrams (e) and (f) 



The two diagrams lead to a vector current 

9A i / k\A ^tt 
-k 

UB\'- 



+ ^{B\ T ^\A)u f T^{q,K)u i . (lor,, 



Here, -P^(g, fc^) is defined in Eq. f )46|) . Ty 7r (g, fc„.) is defined in Eq. fH2|) . and u ^ 1. Finally, 
we can combine diagrams (c), (d), (e), and (f) to get vector current conservation. 
For the axial-vector current, diagram (e) does not contribute, and we find 



+ ^ (S| y |A) 4k, «, + ^—-j ) «, (107) 



+ t^ 4fc Ml-^ + ^J^' (108) 

Here, r^ 7r (g, fc^) is given in Eq. (|33|) and leads to a z/ ^ 1 contribution. The contributions 
due to kv, ^7,-, and K\ are at z/ = 2. It is easy to check that PCAC holds. 

For the baryon current, diagrams (e) and (f) do not contribute at order v = 1: ( Jg) n = 0. 

For the NC photon production matrix element we find 

+ {t) A b 2M ^!l; ) ' , "'W(» ■ (no) 

It is straightforward to see that PCAC is satisfied. Here v = 3; for z/ < 3, there are no contact 
vertices contributing to the NC photon production channel. By power counting, we expect 
that at low energy, these terms can be neglected compared to the v — 1 terms. However, 
as claimed in Ref. (67[, these contact vertices have possible high-energy extrapolations due 
to the anomalous decays of the u> and p. According to Ref. [67J, these terms may play an 
important role in coherent photon production. However, we must realize that the constants 
ci and ei can only be fixed by experiment, and it is not clear that only anomalous meson 
decay will contribute to these operators at low energy. As shown in the lagrangian, we can 
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construct meson dominance by coupling mesons instead of photons to the vertex. Moreover, 
these terms are also the same as operators induced by the off-shell parameters in the A 
lagrangian. 

IV. RESULTS 

In this section, after introducing the kinematics, we will discuss our results for CC and 
NC pion production, and also NC photon production, and compare them with available data 
whenever possible. Aiming at the excessive events in the MiniBooNE experiment, we will 
focus on the scattering of and off nucleons with E v ^ ^ 0.5 GeV. 

A. Kinematics 



Figure [2] shows the configuration in the isobaric frame, i.e., the cm frame of the final 
nucleon and pion. The momenta are measured in this frame, except those labeled as p L , 
which denotes a momentum measured in the Lab frame with the initial nucleon stationary. 
Detailed analysis of the kinematics is given in Appendix [Gj The expression for the total 

2 

cross section for this process is (|M| is the averaged total interaction amplitude squared.) 




FIG. 2: The configuration in the isobaric frame. 



a = 



1 



\M\ 2 



(27T)*tfW(5>) 



d z vlj d% L d 3 p n L f 



frc)*2Ek (27T)32£tf (2tt)3 2 ^ 




2 




(111) 



Based on the equations in Appendix [G], we can make the following estimates: 
For CC pion production: 
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• When Efj = 0.4 GeV, (M nn ) max ~ 1.17 GeV, Q 2 ^ « 0.2 GeV 2 . 

• When El = 0.5 GeV, (M nn ) max ~ 1.24 GeV, « 0.3 GeV 2 . 

We can see that above -E 1 ^ = 0.4 GeV, the interaction begins to be dominated by the A 
resonance. However, when El = 0.75 GeV, (M nn ) max ~ 1.4 GeV, and higher resonances, 
for example Pn(1440), may play a significant role. The exception is that for + p — > 
H~ + p + vr + , only 1 = 3/2 can contribute, and the next resonance in this channel is the 
A(1600), which is accessible only when El ^ 1.8 GeV. 
For NC pion production: 

. When El = 0.3 GeV, (M nn ) max ~ 1.2 GeV, Q 2 max « 0.1 GeV 2 . 

• When El = 0.5 GeV, (M wn ) max ~ 1.35 GeV, Q 2 max « 0.3 GeV 2 . 

Here, we can see that above El = 0.3 GeV, the interaction begins to be dominated by the 
A. However, when El = 0.6 GeV, (M^max ~ 1.4 GeV, and higher resonances may play a 
significant role. 

For NC photon production (E% ^ 0.2 GeV): 

• When El = 0.3 GeV, (M jn ) max « 1.2 GeV, Q 2 max « 0.1 GeV 2 . 

• When El = 0.5 GeV, (M yn ) max « 1.35 GeV, Q 2 ^ w 0.3 GeV 2 . 

Here, we expect the A to dominate when El ^ 0.3 GeV. But, similar to the case of NC 
pion production, higher resonances may need to be considered when El ^ 0.6 GeV. 

From the analysis outlined above, we can expect our EFT to be valid at El ^0.5 GeV, 
since only the A resonance can be excited, and Q 2 ^ 0.3 GeV 2 , so that meson dominance 



works for various currents' form factors [20f| . To go beyond this energy regime when we show 
our results, we will require M nn 1.4 GeV and will use standard phenomenological form 
factors that work when Q 2 ^ 0.3 GeV 2 . 



B. CC pion production 

In this section, we will compare CC pion neutrinoproduction results with ANL [HI and 
BNL [35j measurements. In both experiments, the targets are hydrogen and deuterium. (All 
the other experiments use much heavier nuclear targets in (anti) neutrino scattering, and to 
explain this, we must examine many-body effects.) The beam is u^, the average energy of 
which is 1 GeV and 1.6 GeV for ANL and BNL, respectively. In the ANL data, there is a 
cut on the invariant mass of the pion and final nucleon system: M vn ^ 1.4 GeV. In the 
BNL data, there is no such cut. Based on the phase-space analysis discussed above, this cut 
clearly reduces the number of events when E v is above 0.5 ~ 0.6 GeV. Since the data stretch 
above this limit, in the first three figures: EJHJ and|5l we show our conventional form factor 
( £ cff ') calculations with the M^ n constraint. That is, for F md , G md , ca, and d& we substitute 
the conventional form factors used in the literature j32[. Then we apply our lagrangian to 
the meson dominance form factor ('mdff ') calculations. As we have already concluded that 
this 'mdff' approach is inadequate above E v = 0.5 GeV, in the following figures: El [TJ El El 
E3 andHH we show the 'mdff' results with E v < 0.5 ~ 0.6 GeV, for which M nn sC 1.4 GeV 
holds automatically. Since we believe the EFT is applicable in this low-energy regime, in 
these figures, we show results including Feynman diagrams up to order v = 1 and v = 2. 
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FIG. 3: Total cross section for v^+p — > \i~ + p + n + . 'Only A' indicates that only diagrams with 
A (both s and u channels) are included. 'Up to v = V includes all the diagrams at leading order. 
The code 'cff' indicates that the calculations are done with conventional form factors, while 'mdff' 
indicates the calculations are based on the EFT lagrangian with meson dominance. In the ANL 
data, M m ^ 1.4 GeV is applied, while no such cut is applied in the BNL data. For all calculations, 
M nn ^ 1.4 GeV is applied. 



In Fig. [31 we show the data and calculations for + p — > p~ + p + tt + . The ANL data 
is systemically smaller than the BNL data, due to enforcing the M nn constraint at ANL. 
As mentioned above, we make use of the conventional form factors and include in the 'cff 
only A' calculation the Feynman diagrams with the A in both s and u channels and in the 
'cff up to v — V all the Feynman diagrams up to leading order. These two calculations are 



quite similar to those done in Ref. [30[. Indeed, our results are consistent with theirs for 
the conventional value of (In Ref. (30[, only the s channel contribution is included in 
the calculation with 'only A'.) Next, we apply our lagrangian in the 'mdff' calculations, in 
which form factors are realized in terms of meson dominance. In Fig. [3j we show both the 
result with only A diagrams and the result with all the leading-order diagrams in the 'mdff' 
calculations, so that we can compare the 'mdff' approach with the 'cff' approach. 

First, we can see that both 'cff' and 'mdff' with only A diagrams are consistent with the 
data at E v ^ 0.5 GeV. Introducing other diagrams up to order v = 1 is still allowed by the 
data at low energy, although they indeed increase the cross section noticeably. Second, the 
two approaches with the same diagrams begin to differ from each other beyond E v = 0.5 GeV, 
which is also consistent with the analysis of phase space and the discussion of the validity 
of meson dominance. In Ref. [30], a reduced C^(0) is introduced, primarily to reduce the 
calculated cross sections above E v — 1 GeV. However, since we are only concerned with the 
E v ^ 0.5 GeV region, in which we see satisfactory agreement between our calculations and 
the data, we will stick to the Cg (0) fitted from the A free width in the framework of our 
EFT lagrangian. Furthermore, in the original spectrum-averaged da/dQ 2 data of ANL 34| . 
the contributions from E v ^ 0.5 GeV neutrinos are excluded, so comparing calculations with 
data at low energy is not feasible at this stage, and we will not show our da/dQ 2 here. 

In Figs. S] and El we show the data and calculations for + n — > p~ + n + tt + and + 
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n — y /i~ +J9 + 7T . We can see that the situations in these two processes are quite similar to 
the one in Fig. 0J the results of the 'cff' and 'mdff' approaches are consistent with the data at 
low energy. Again the differences between the two approaches with the same diagrams begin 
to show up when the neutrino energy goes beyond 0.5 GeV. Although the pion production is 
still dominated by the A, other diagrams introduce significant contributions, which violate 
the naive estimate of the ratio of the three channels' cross sections based on isospin symmetry 
and A dominance. 

In Figs. El El El Ell and dTJ we begin to investigate the convergence of our calculations 
in different channels in neutrino and antineutrino scattering. We show the 'mdff' calculations 
based on our EFT lagrangian up to different orders. We see that the power counting makes 
sense systematically in different channels: including N and contact terms up to v = 1 
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FIG. 7: Total cross section for i/„ + n — > fi + n + ir + . The curves are defined as in Fig. [H 

changes the 'only A' calculation non-negligibly. (Far away from resonance, we see that the A 
contribution is not dominant compared to other diagrams, and it begins to dominate around 
0.4 GeV. This is consistent with the power counting discussed in Sec. IIIIBI) . However, the 
v = 2 terms do not change the 'up to v = V results significantly. This partially justifies 
the use of meson dominance, which automatically includes higher-order terms. All the 
calculations of neutrino scattering are consistent with the limited data from ANL. We can 
see that the cross section for antineutrino scattering is generally smaller than that of neutrino 
scattering, due to the relative sign chosen between V tfl and A ltl in the Feynman diagrams 
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with A. The signs between V 1 ^ and A 1 ^ in other diagrams is well defined in our lagrangian. 
However, the relative sign between currents due to the A and other diagrams is not well 
constrained by the available data in our framework, as indicated in Ref. [30(]. Here we rely 
on the sign of Ka fitted in pion-nucleon scattering |68| to set the sign between the A current 
and the background contribution, as shown in Eqs. (171]) to (1501) . 

°- 45 I 1 tjptov=2 l n 1 ] 
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FIG. 8: Total cross section for + n — > fj,~ + p + 7r°. The curves are defined as in Fig. [6l 
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FIG. 9: Total cross section for Vn + p — > fi + + p + ir~ . The curves are defined as in Fig. El 



C. NC pion production 

In this section, we discuss the results for NC pion production in (anti) neutrino scattering. 
In Figs. [12] and [131 the results in the 'mdff' approach including diagrams of different orders 
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are shown for neutrino scattering, while the results for antineutrino scattering are shown in 
Figs. [14] and [151 The channels are explained in each plot. 

Since all of the available data for NC pion production are spectrum averaged, and neutri- 
nos with E u ^ 0.5 GeV have small weight in such spectrum integrated analyses, we will not 
compare our results with data. In other words, current data does not put strong constraints 
on the NC pion production in this energy regime. 

Nevertheless, we clearly see the convergence of our calculations; introducing the v — 2 
terms does not change the total cross section significantly. However, we also see the violation 
of isospin symmetry in the 'up to v = V and 'up to v = 2' calculations in each plot, if we 
compare each pair of channels in every plot, namely, Figs. [T3J [HJ and (TBI In principle, 
if there is no baryon current contribution in NC production, we should see that the two 
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FIG. 13: Total cross section for NC ir^ production due to neutrino scattering. The curves are 
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channels yield the same results in each plot. For example, isospin symmetry implies 

(p,7c°\V°^A^\p) = (n,7c \V°^A 0lM \n) , (112) 

(p,Ar B \p) = -{n,n°\J» • (H3) 

So with 'only A', we will not see the difference between the two cross sections, since the 
(isoscalar) baryon current cannot induce transitions from N to A. After introducing back- 
ground terms, which contain contributions from the baryon current, we would expect the 
results for the two processes to be different, as confirmed in Fig. [12j This analysis is appli- 
cable to the other plots. 
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D. NC photon production 

In this section we focus on NC photon production. Besides NC 7r° production, this process 
is another important background in neutrino experiments. As we are ultimately concerned 
with the excessive events in the MiniBooNE experiment [l[, we focus on E u ^ 0.5 GeV, 
as mentioned in the beginning of this paper. One important difference between NC photon 
production and CC and NC pion production, is that all of the v — 2 terms do not contribute 
in this process. Therefore, we include the two v = 3 terms in NC photon production, 
namely, the e\ and c\ couplings in Eq. (IllOp . As mentioned in Sec. Ill B 1\ there are many 
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other interaction terms at the same order as the e\ and c\ terms, but these are the only 
two contributing in NC photon production. Moreover, these two couplings are singled out in 
Ref. |6_7| as the low-energy manifestations of anomalous p and u decay and are believed to 
give important contributions in coherent photon production from nuclei. Here we investigate 
the consequences of these two couplings. We emphasize that from the EFT perspective, the 
only way to determine these two couplings is by comparing the final theoretical result with 
data, rather than by calculating them from anomalous decay, which is not necessarily the 
only higher-energy physics contributing to these two operators. For example, as we discussed 
the off-shell couplings before, an off-shell coupling between N, tt, and A can introduce the 
same matrix element as that induced by these two contact terms. Changing the off-shell 
couplings would also change these two contact terms to make the theory independent of 
the choice of off-shell couplings. Nevertheless, to perform concrete calculations with these 
two terms without precise information on the coupling strengths, we use the values from 
Ref. (H3 in Figs. EH and ED 

We can see the convergence of our calculations. The two couplings introduced in the 'up 
to v — 3' calculations increase the total cross section in both channels for both neutrino and 
antineutrino scattering, although the change is quite small. This constructive behavior is 
consistent with the results in Ref. [67]. However, as is easily seen, the contributions of this 
process are negligible in scattering off a single nucleon. 
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FIG. 16: Total cross section for NC photon production due to neutrino scattering. 'Only A' 
indicates that only diagrams with A (both s and u channels) are included. 'Up to v = V includes 
all the diagrams at leading order. 'Up to v = 3' includes higher-order diagrams. The v = 2 terms 
are zero in these channels. The next-to-leading order in this channel is v = 3, whose couplings are 
from Ref. 67]. 



Naive power counting, however, does not give an accurate comparison between the A 
contributions and the N contributions at low energy. One reason is that the neutron does 
not have an electric charge at low energy, so its current should appear at higher order than 
the naive estimate. The second reason is that for the proton, due to the cancelation between 
the baryon current and the vector current, the neutral current is mainly composed of the 
axial- vector current, which reduces the strength of the neutral current. Because of these two 
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FIG. 17: Total cross section for NC photon productions due to antineutrino scattering. The curves 
are defined as in Fig. [TBI and the channels are also indicated. 

factors, the contributions of A and N are at the same scale far away from resonance, but 
near the resonance, the A dominates. 

V. SUMMARY 

Weak pion and photon production from nucleons and nuclei produce important back- 
grounds in neutrino-oscillation experiments and must therefore be understood quantita- 
tively. In this work, we studied pion and photon neutrinoproduction in a Lorentz-covariant, 
chirally invariant, meson-baryon EFT. For neutrino energies E\f° < 1 GeV, the resonant 
behavior of the A is important. We therefore included the A degrees of freedom explicitly 
in our EFT lagrangian, in a manner that is consistent with both Lorentz covariance and 
chiral symmetry. 

It is well known that in a lagrangian with a finite number of interaction terms, including 
the A as a Rarita-Schwinger field leads to inconsistencies for strong couplings, strong fields, 
or large field variations. In a modern EFT with an infinite number of interaction terms, 
however, these pathologies can be removed, if we work at low energies with weak boson fields. 
This is because the problematic terms in the lagrangian produce local contact interactions 
that can be absorbed into other contact terms in the EFT lagrangian. Ambiguous, so-called 
off-shell couplings have also been shown to be redundant in the modern EFT framework. 
Thus the A resonance can be introduced into our EFT lagrangian in a consistent way. 
Moreover, we studied the structure of the dressed A propagator and found that it has a pole 
only in the spin-3/2 channel, so that we indeed have the correct number of resonant degrees 
of freedom. 

Because of the symmetries built into our lagrangian, the vector currents are conserved 
and the axial-vector currents satisfy PCAC automatically, which is not true in some of the 
other approaches to this problem. Needless to say, a conserved vector current is crucial for 
computing photon production. By using vector and axial-vector transition currents that 
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were calibrated at high energies, we found results for pion production at lower energies that 
are consistent with the (limited) data. This was also true when vertices described by meson 
dominance were used. We also studied the convergence of our power-counting scheme at low 
energies and found that next-to-leading-order tree-level corrections are very small. Finally, 
we computed neutral-current photon production including contact interactions consistent 
with anomalous p and u decays and found that, at least for a nucleon target, the resulting 
cross sections are unmeasurably small. 

We are currently using this QHD/EFT framework to study the electroweak response 
of the nuclear many-body system, so that we can extend our results to pion and photon 
neutrinoproduction from nuclei, which are the true targets in existing neutrino-oscillation 
experiments. 
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Appendix A: isospin indices, T matrices 



Suppose t are the generators of some (ir)reducible representation of 577(2); then it is easy 
to prove that, in matrix form (6 = — e~ mtv ), 



6t5 



->T 
-t , 



(Al) 



where the superscript T denotes transpose. This equation justifies the use of 5 as a metric 
linking the representation and the complex conjugate representation. One easily finds for 
V^ 2 \ DM and^ 1 / 2 ); 



cab 



( l\ 
0-10 
10 
\-l 0/ 



>ab 



/o -l\ 

10 
0-100 

\1 J 



5 ij 
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AB 




JAB 




We turn now to the T matrices. As discussed in Sec. Ill 

rrr\ iA 



2' ^ 



T a lA = (l,~;i,A||;a) . 



(A2) 

(A3) 
(A4) 



(A5) 
(A6) 
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To be more specific: 
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(A7) 



a A 



1 



Aa 



It is easy to prove the following relations (here t 1 is a Pauli matrix): 



= 5 i + ze ^ rfc > 



Aa 

(A8) 

(A9) 
(A10) 
(All) 



Here P{ is a projection operator that projects £g> 'H^ 1 -' onto T^a). 

A few words about e"*- fc are in order here. We have the following transformations of pion 
fields: 



TT = UtTC 



here, i = +1, 0, — 1 ; I = x,y,z ; 
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(A12) 



(A13) 



and hence we have 



_ _ _X„X I —11— .1/ I —Z„Z 

= -TT + V- 1 - TT- 1 ^ 1 + 7l Z 7T Z 
= 7T +1 7T +1 + 7T~ 7T_i + 7r°7To . 



(A14) 



Meanwhile, under such transformations, 



rijk 



fijk 



V^ T U j rU 1 le IJK 



i \ ijk 



j-uj-u k c — det(u))e 

—i, if ijk = +1, 0, — 1 ; 
-i6 v , if ijk = P(+l, 0, -1) . 



(A15) 



Appendix B: C, P, and T symmetry realized in QCD 



The symmetries C, P, and T are preserved in QCD. Based on the lagrangian in Eq. (J5J) 
we find the corresponding transformation rules shown in Table [B Inside the table, Vjj = 
diag(l, -1, -1, -1)^ and Tf = diag(-l, 1, 1, 1)^. 
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TABLE I: Transformations of background fields under C, P, and T operations. The transformed 
spacetime arguments are not shown here. 



Appendix C: C, P, and T symmetry realized in QHD 



The C, P, and T transformation rules are summarized in Table [III A plus sign means 
normal, while a minus sign means abnormal, i.e., an extra minus sign exists in the transfor- 
mation. The convention for Dirac matrices sandwiched by nucleon and/or A fields are 



CNTNC- 1 
C(ATN + NTA)C- 1 
Ci(ATN -NTA)^ 1 



N T T T N T 



-A T T T N T 



normal ; 
abnormal 



N T T T A T 



-A T T T N T + N T T T A T 



/A' V N — iN r A 
-tA T T T N T + iN T T T A T 



normal ; 
abnormal . 

normal ; 
abnormal 



(CI) 
(C2) 
(C3) 



Here, in Eqs. (101 j) . (102|) . and (103|) . the extra minus sign arises because the fermion fields 
anticommute. The factor of i in Eq. (1031) is due to the requirement of hermiticity of the 



lagrangian. To make the analysis easier for A YN + C.C., we can just attribute a minus sign 
to an i under the C transformation. Whenever an i exists, the lagrangian takes the form 
i(ATN — NT A). When no i exists, the lagrangian will be like A TN + NT A. 

For P and T transformations, the conventions are the same for A" and A fields, except 
for an extra minus sign in the parity assignment for each A field j69[, so we list only the A" 
case (Vjf and TJ 1 can be found in Appendix [B]) : 



PNT^NP- 1 



TNT^NT- 



nv;t u n , 

-NV^T U N 



NT U U T U N 



normal ; 
abnormal 

normal ; 



-N 77 r„ N , abnormal . 



(C4) 
(C5) 



It is easy to generalize these results to T^, etc. 

Suppose an isovector object is denoted as M = O^t 1 , then the conventions are explained 
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TABLE II: Transformation properties of objects under C, P, and T. 



below: 

cr 1 = ' normal ; (C6) 

[ — , abnormal . 

P0 p-i = J ' normal 5 ( C7 ) 

-V^O u , abnormal . 

TO^ = { T »°: n 7 mali , (C8) 
— / ' U v , abnormal . 

The same convention applies to the isovector (pseudo)tensors. For isovector (pseudo)scalars, 
the V and T should be changed to 1. For the C transformation, T means transposing both 
isospin and Dirac matrices in the definition of O, if necessary. 



Appendix D: Form factors for currents 

Here we use matrix elements of the various currents to define the form factors produced 
by the EFT lagrangian j^J. Note that is defined as the incoming momentum transfer at 
the vertex; in terms of initial and final nucleon momenta, = p^f — Pni- When a pion is 
emitted, p^ + q^ = p^ f + k%. 

To derive these expressions, it is necessary to expand terms in the lagrangian to leading 
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order in pion and external fields. Useful results are given below: 
Tr&[U,&>rf]) « 2^ fc ^^, 

^ J TV J 71 

i i 

Tr(l-{U « -2i ' 



2 J/ /. ' 



? 2^ 2 2 ' 

? 2 ? 2 ^ 2 ' 



~ ijk a r « r ijfc r fc 

~ 2/2 e 7r .? c V rfe "2 ~ Vifl ~2~ e y~~2~ a<At ' 

+ background interference terms, 

/l/b/ + /r^ = 2<9 [At v H - 2?[v M , v„] - 2i[a M , a„] , 
fLvv - fR»v = -2d [fM a u] + 2i[v M , a„] + 2i[a M , v„] , 

A Jtt 

~ 29[ M v^] + 2e lJ ' fc -p — ^ajj,] + background interference, 

/•7T 2 

"TrO " Li{iv fllifiv) ^ ^ "T - 7~k{fLi^u fuifiu) 
A Jtt 

~ — 2<9[ M aj,j — 2e lJ ' fc -r- — <9[ M v^] + background interference. 

J7T 2 

We now proceed to determine the matrix elements. 
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(N,B\V;\N,A) = 
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= -j 1 (B\ — \A} u f T Anfl (q,k 7T )u i 



(D23) 



Now we consider (N, B\A l \N, A) and (N, B; ir, j\V*\N, A). In the chiral limit, we find 



T 



(N,B\A;\N,A) = -(B\-\A}u f Yl 5 



U; 



9 A [ 9^u ^ 



1 (5 9^ - %<lv) 



M 2 



9 q Qjiv - 



(D24) 



(A^;7r,jA|^iA^> 



9a9^v ~ jp [q ■ (q - K)9ixv - (q - K)nQu] 



- 2c ai g a 



q-(q- K)9fiu - (q - K)uQ. v 



(q - k w ) 2 - m\ 



ai 



(D25) 



Now suppose there is only one manifestly chiral-symmetry-breaking term, i.e., the mass term 
for pions; then the pion-pole contribution associated with the 9a coupling in (TV, B\A^\N, A) 
will become 9a[9hv — q^qu/{q 2 — while the other parts in (N, B\A l \N, A), as well as 

the whole {N, B; ir,j\V*\N, A), will remain unchanged. However, we must realize that there 
are other possible chiral-symmetry-breaking terms contributing to (N, B\A l \N, A). For 
example, ijnl/M) Nrf(U - U^)N will contribute to (N, B\A^\N, A) as 

_ 2ml q» 4i b iB \ Ti \ A \ 

M 2 q 2 — m 2 1 2 1 ; ' 
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To simplify the fitting procedures, we will use the following form factors: 



(N,B\AUN,A) = -GT{q A ){B\-\A)u } [ gilv 



2 \ q 2 — m 2 



(D26) 



,- I 



(N,B-ir,j,K\V;\N,A) = ^ (B\-\A) u f Yj 5 Ui 



j ^ Q ■ (g - K)g^ - (g - K)»gv 

(q ~ K) 2 



+ SGT((q-K) 2 ) 



The required definitions can be found in Eqs. ( 1431) and ()44l) . 
Finally, we calculate the pion form factor (tt, k\V^\ir,j): 



(n,k,k n \V*\7r,j,k w -q) = ie\(2k w - g) M 

+ 2i 



9pTTiT ij q i 



/ ~> „2 _ ,,,2 ■ q 2 ^^) 



q 2 — > m 2 in numerator 



# 7 g^ - m p 



+ 2^ ii 



I, o _ ,.,2 ' ^(fc ? 2 ^m) 



g 7 q* — m 



pion on shell = ie 1 ^ (2^ — g) M I 1 



9pTTir q 

g 1 q 2 - m 2 p 



ie\ (2K - g)^^ (g 2 ) 



(D27) 



(D28) 



(D29) 



(n,k,K\V!\n,j,K-q) = ie 



(2^ - g), + 2^(9") ( ^ - ^ g„ 



Here, M^V) 



ie lJ k P Vfl (q,k n ) . 
F™ d {q 2 ) - F™ d (0) . 



(D30) 
(D31) 
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Appendix E: Free A propagator and self-energy insertion 



Normally, the free, spin-3/2 field's propagator can be decomposed as |45l. l7d. 171 



sf v ip) 



i 



( I 



— m + ie 
2 



V3> 



P 



m 



12 



p(j)p 

Jo 



21 



3m 2 
1 
3 



(^+m)P 2 ( J^ 



1 

3p 



, 7 L T 



3p 



1 



21 



-*22 



V% 2 

1 

TP 
i 



It is easy to prove the following relations: 

(p CO r(p (i) 



p® + Pip + p$ ) 



Pi 



(I) 



11 



22 
22 



p(D , p' 



P 



(I) 

11 7 
p(|) 

-^22 ) 



r iw y 

pjJV 
1 , 

pCf-L) 

0, 
0, 

. 








(El) 
(E2) 

(E3) 
(E4) 
(E5) 
(E6) 

(E7) 
(E8) 
(E9) 
(E10) 
(Ell) 
(E12) 

(E13) 
(E14) 
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With the preceding properties of the projection operators, we easily find 

-1 „^ 1 „(D 1 „rii _(i) 



s°Ap) 



p(§) 



p(f) 



^ — m + ie 

-1 

— m + ie 

1 



p(§) 



v/3 



m 



^12 



\/3' 



p( 2 ' _ p 2 

^21 i -<22 



2 (^ + m)P 2 ( | } 



3m : 



p(f) 



v3m 

,0(2 



p(|) 



\/3' 



pW I 0(2) 

-'21 "T 



22 



3m 2 



+ m)P 2 ( | ) 



P^) (E15) 



(E16) 



Following the analysis in Ref. [45[, the self-energy of the A can be defined as E 



S A j ^ + It follows that in tfE^, the indices can only have structures like the prod- 

Then we find quite an interesting property of S (the fiu indices are 



UCtS (jv,P»)(jv,Pv) 

suppressed) : 



= (p(D + pd ± ))(S A ( 7 + 5E)(P<f> + pd^) 
= P®E A gP® + pd^SpCf- 1 -) 

+ P(§)(£ A £ + 5S)pd ± ) + p(i ± )(E A 5 + 5S)p(i) 

= p(Ds A p(f) + pd^spd^ 

= E (|) + E (|i) . 



(E17) 



In the proof, we make use of 



p(f) 



and 



p(§±) E A 



0, since the only possible 



spin structures of S A are 1, p 1 , and 7 5 (parity violation), which commute with the two 
projection operators. This implies P^E^P^- 1 ) = and P^-^E^P^) = 0. We also 
make use of Eqs. (|E8]) and (jE9|l . so we get pd^EPd^) = and pd^EPd ) = 0. 



Appendix F: Construction of a Delta interaction term 



We show an example of constructing a term in the lagrangian with interactions between 
pions and the A. Consider 

^a;{^,^ 7 5 }a 6 ,. 

It can be shown that 

= A a » ^7 5 A bM + A; [-Ya^-Ya^ + Y ^tV]'** • 

(Fl) 

According to the argument in Sec. Ill C 3\ the original coupling and A aM ^ 7 s A bfl are equiv- 
alent in low-energy effective theory. We use the second form in our lagrangian. 
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Appendix G: kinematics 



Following a standard calculation, we find the total cross section: 



a = 



Mpfi-Pflil 1 ) YT 1 J (2^2^ (2tt)32^ (2tt)32^ 



/ 



/ 
/ 



Mpfi-P L m\ 1 J l " Pm ^ ^(2^)32^(2^)32^(2^)32^ 
4bi " P&| 1 J [q Pm Pnf PJ (2tt)32K/ (2tt)32S§ (2tt)32^ 



2 

M 1 






32M„ (2vr 


) 5 K + £„/ 





dn n dEfjdnf f . (Gl) 



It is quite complicated to calculate the boundary of phase space in terms of the integration 
variables in the preceding equations. Later, we will work out the boundary of phase space 
in terms of the invariant variables Q 2 and M nn in the cm frame of the whole system, so we 
would like to have the following: 

= _ M 2 + 2E L {E L _ CQS 9 L f) j (G2) 

Ml = (q L + p L ni ) 2 = -Q 2 + M 2 n + 2M n (Et - E^ f ) , (G3) 
from which it follows 

dQ 2 dM 2 n = AM n Et \tff\dEfidcoseff . (G4) 

By using the invariance of the cross section with respect to rotations around the incoming 
lepton direction, we have J dilf f = J dcosd^ 2n, and thus 

- 1 Sik^wk^'^ 2 ■ (G5) 

In the isobaric frame, there is no preference in the direction of the outgoing pion. Thus 
the boundary of Q n is the whole solid angle in the isobaric frame. Now let's work out the 
boundary of phase space in the cm frame. We have 

M\ = v\= ip L m + Ph? = (M n + Et) 2 - {Eli) 2 = M 2 n + 2MX , (G6) 
Ml n = (p, + p r ,f) 2 = ( P C A - P? f f = M 2 A + M? f - 2M A E? f . (G7) 

Here, Ef^ is the final lepton's energy in the cm frame. From now on, all the quantities in 
the cm will be labeled in this way. So, for given E^, i.e., Ma, we can see that 

M„ I M, < M wn <: M A - M lf . (G8) 
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By using Eq. (1G7j) . we find 



Then, for given and M™ (or using Q 2 = -Mf f + 2E%Ef f -2E%\p lf \ c cos 6^ (where 

in 1~ i-i /~\ . "v i r\* I /"i KArnmnn r Kin nnr rri - 



6*£ is the angle between the outgoing lepton's direction and the incoming lepton's direction 



in the cm frame, and Ef] = {M\ — M 2 ) /2Ma is the initial lepton's energy in the cm frame), 
we finally arrive at 

[Q\Ef f )] min = -Mf } + h lf , (Gil) 

[Q 2 (Eg)] max = -Mj + 2Eg (JS§ + yj (E°)* - M? f ) . (G12) 

These equations give a description of the phase-space boundary in terms of the invariants 
M nn and Q 2 . 
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